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We study the weak-tunnelling limit for a system of cold """K atoms trapped in a one-dimensional optical lattice 
close to an 5-wave Feshbach resonance. We calculate the local spectrum for three atoms at one site of the lattice 
within a two-channel model. Our results indicate that, for this one-dimensional system, one- and two-channel 
models will differ close to the Feshbach resonance, although the two theories would converge in the limit of 
strong Feshbach coupling. We also find level crossings in the low-energy spectrum of a single well with three 
atoms that may lead to quantum phase transition for an optical lattice of many wells. We discuss the stability of 
the system to a phase with non-uniform density. 
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I. INTRODUCTION 

Quantum simulation of strongly correlated models of con- 
densed matter systems has become one of the most exciting 
goals in ultra-cold atom experiments, thanks to the spectac- 
ular progress in trapping cold fermions in optical lattices at 
degeneracy temperatures This new system goes beyond 
traditional experiments in solids, because it is readily possible 
to tune independently the hopping amplitude and interaction 
strength. In particular, by changing the external magnetic field 
that tunes the relative Zeeman energy between the scattering 
state and a closed (bound state) channel, one can reach a par- 
ticular Feshbach resonance between two hyperfine states of, 
for example, "^''K where the scattering length between these 
states becomes very large and positive [2|. This then allows 
one to tune the on-site two-body interaction to be strongly re- 
pulsive. With this, the fermionic Mott insulator, where inter- 
action driven locking of fermions into a crystalline structure, 
has been demonstrated in two recent experiments 13] |4l . Cur- 
rently, there are great experimental interests to reach ordered 
phases in these systems, in order to understand the phase di- 
agram of the two-dimensional single-band Hubbard model, a 
model that has been hypothesized to be the minimal model for 
the high temperature superconductors 0. 

Indeed, the unprecedented experimental possibility of 
studying extremely strong repulsion has spuiTed on much new 
theoretical effort |jT]. In particular, it was recognised that, 
close to a Feshbach resonance, the strongly enhanced two- 
body interaction leads to states with occupation of multiple 
Bloch bands of the optical lattice, both in experimental ||6| and 
theoretical studies |7, 8|. These multi-band systems 191 ITOl 
may lead to richer phases, by analogy with multi-orbital sys- 
tems in solid state materials. 

While the ultra-cold atoms parameters of the one-band 
Hubbard model have been derived right from the beginning 
0T|, the parameters for the multi-band generalisation have 
only more recently received some attention [7, J] [12] [T3] [14] 
[T5l . Motivated by the experiments of Kohl etal. ||6l where the 
average filling per lattice site is about two fermions or slightly 
more, we study here the energy spectra when two or three 
fermions interact near a Feshbach resonance within a site of 
the deep optical lattice, using a two-channel description of the 



Feshbach resonance. This is then the first step towards de- 
riving from microscopic physics, the effective model and its 
parameters for a system of strongly interacting fermions in an 
optical lattice. 

In this paper, we consider a system of cold fermionic "^"K 
atoms in the two lowest internal states trapped in a deep op- 
tical lattice close to an inter-species i-wave Feshbach reso- 
nance. We assume a lattice with average filling of two atoms 
per site with an equal mix of the two internal states. For 
simplicity, we consider only a one-dimensional (ID) system. 
In experiments with cold atoms in optical lattices, a set of 
quasi-lD tubes can be obtained from a 3D optical lattice, 
V(r) - Voi- sin^ k^x + Vq,. sin^ fcyj + Vo,- sin^ k,z, by adjust- 
ing the laser amplitudes so that Vqz ^ Vqx = Voy lfT6l [TTl . 
We concentrate on the weak tunnelling limit (deep lattice), 
in which experiments with lattice fermions across a Feshbach 
resonance are usually performed |6 |, and approximate the op- 
tical lattice potential at each lattice site by a harmonic well, 
= 2i=x,y,z/''^^/2 with - 2kfVoi- The local spec- 

trum and eigenstates for two and three atoms at a lattice site 
provide then a microscopic basis for the model of such lat- 
tices found by Ho ifTOl where a weak tunnelling is introduced 
perturbatively. 

The spectrum for two atoms in a harmonic well close to 
a s-wave Feshbach resonance was previously studied in e.g. 
Q within the one-channel model for the harmonic confine- 
ment in one, two and three dimensions and in f8^| within the 
two-channel model for a three dimensional harmonic trap. 
However, the spectrum for three atoms in a well was previ- 
ously studied only within the one-channel model (e.g. for 3D 
harmonic confinement l[T2l , quasi- ID confinement |[T3]| and 
fermions at unitarity |T4l [TSl ). Here, we calculate the spec- 
trum within the two-channel model and examine the ground 
state for our system across the resonance. We find the level 
crossings in the low energy spectrum that may lead to quan- 
tum phase transition when multiple wells form an optical lat- 
tice. Similar features were found in the spectrum for a 3D 
system in [121 within the one-channel model. However, we 
find that there are quantitative differences between our two- 
channel and their one-channel results at and near the Feshbach 
resonance. 
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II. TWO-CHANNEL MODEL FOR THREE ATOMS IN A 
WELL 

When Voz ^ Vqa = Voy, the system can be described by an 
effective ID Hamiltonian. The effective ID parameters can be 
related to the physical 3D parameters in the limit — > as 
presented in 1 18 1. We concentrate on the weak tunnelling limit 
and calculate the local spectrum when there are three atoms in 
two different internal states in the well. 

The "^"K atoms are assumed to be trapped in the two 
lowest hyperfine energy states \F-i,mf^} = |9/2,-9/2) and 
\F2,mf^} = \9 12,-1 12). The interatomic interaction poten- 
tial, Uifx - i*2) = Uciri - ri) + UAri - r2)^i ■ Si, de- 
pends only on electronic spins of atoms and preserves the to- 
tal angular momentum projection of the two colliding atoms, 
nif - nif^ + nif^. Then, for the two ^"K atoms in the two 
lowest internal states that constitute the open entrance chan- 
nel, \ai,a2) = 19/2,-9/2; 9/2, -7/2), the only closed chan- 
nel coupled to this state is |9/2, -9/2; 7/2, -7/2). For sim- 
plicity, we will use here the notation |9/2, -9/2) s | f), 
|9/2, -7/2) = I i) and |7/2, -7/2) = |o). In this notation, the 
Feshbach physics relies on the scattering between the states: 



dinates: 



T°) 



(1) 



where the final closed-channel state is a tightly-bound | f o) 
dimer We are interested in a regime near a Feshbach reso- 
nance so that the Feshbach-induced interaction is strong. We 
will therefore neglect for simplicity the direct i-wave scatter- 
ing of the open-channel states. 

In this section, we will focus on the case where there are 
two atoms of the species t and one atom of the species J, in 
the well. (We will comment on the case of one t and two i 
atoms at the end of the section.) We will calculate how the 
open-channel levels | T.T.i) ('-e-. atoms 1 and 2 are t and 
atom 3 is X) become coupled together with levels involving a 
pair of atoms in the closed channel | t, °)- 

Consider first the states in the open channel in the absence 
of Feshbach scattering. Since we are ignoring direct .s-wave 
scattering in the open channel, atoms in the open-channel 
states are non-interacting and they will have a Hamiltonian /!°p 
that consists only of the potential of the harmonic trap /i'"p. 

We will work in the centre-of-mass frame of this three-body 
system because the motion of the centre of mass can be fac- 
tored out completely for a single harmonic well. We denote 
the positions of the atoms T.T.-l with = r^z, ri = r2Z and 
'*3 - f^z, respectively. It is conventional to describe the rela- 
tive degrees of freedom of the system using the Jacobi coordi- 
nates: xi - 2 [r3 - (ri + r2)/2] / V3 and X2 - r2 - ri. We find 
that it is more convenient for the evaluation of our formulae 
to use a rotated version of these coordinates: 

X = ^^^^ = ^[2r3 + (V3-l)ri-(V3 + l)r2] 



V2 V6 

Xi + X2 1 

V2 ~ V6 



[2r3-(V3 + l)ri+(V3-l)r2] (2) 



(v^V?) + ^(.^+/) (3) 



2f, 



where yu - m/2 is the reduced mass for the relative motion of 
two atoms. The characteristic length of the oscillations in the 
harmonic trap is d,- - (S/yUw)'^^ in the centre-of-mass frame. 
The eigenstates of this Hamiltonian must have wavefunctions 
which are antisymmetric under the exchange of two fermions 
in the same internal state (t), i-e. under the exchange ri <-» r2 
which corresponds to x y. Therefore, the eigenstates are 

\ui,(x, dr)ui(y, dr) — Mi.(y, dr)ui(x, dr)] 

\k, I) = ^ *^ '^^ — ' ''^ I T, T, i) (4) 

V2 

for non-negative integers k + I and 

uJs, d) — T— — — ■ H, 

with Hp being the Hermite polynomials of order p - 0,1, 

These states have energy (k + I + l)fi(L). Note that Pauli ex- 
clusion is enforced in that the wavefunction vanishes unless 
ki^l. 

Having found basis states for the open channel, we can do 
the same for the closed channel with atoms 1 and 2 in the t 
state and atom 3 in the o state. They will have an energy rela- 
tive to the open channel of 6b which is tunable by an external 
magnetic field B. Moreover, there is an additional attractive 
interaction V^' in the closed channel that creates the | t °) 
dimer essential for Feshbach physics. Therefore, the closed- 
channel Hamiltonian h'^^ is given by: 



■0 



(5) 



(6) 



where the coordinates r-t - - ri 2 represent distances be- 
tween atom 3 (in o state in the closed channel) and each of the 
other two t atoms. 

To find the eigenstates of this closed-channel Hamiltonian 
h^^, it is useful to identify parts of h^^ which provide an at- 
tractive interaction between one pair of t ° atoms to form a 
dimer while leaving the other t atom as a spectator. Let /j';' 
provide an attraction between atoms 1 and 3, leaving atom 2 
as a spectator and /I'i' provide an attraction between atoms 2 
and 3, leaving atom 1 as a spectator. 

These truncated dimer-spectator Hamiltonians are: 

ht = -| (Vi + Vi) + ^ [si + rl) + y^V.) + 6s (7) 

where s± - [ri_2 + rj, - 2r2,i]/ V3 are proportional to the dis- 
tances between the centre of mass of the dimer and the spec- 
tator 

The eigenstates of these trancated Hamiltonians (|7]| are 
atom-dimer states of the form u„(s-t,dr)xif±,d,n). The first 
part, u„, represents the oscillation in the relative displacement 
between the spectator and the dimer with energy (n + 1 /2)ha>+ 
6b- The second part, xi^±^d.,n), is the wavefunction for the 
bound state formed due to the attractive interaction V^': 



The potential of the harmonic trap is separable in these coor- 



V, 

2^l 



(8) 
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with energy ei,. We cannot solve the bound-state equation in 
general. We will simply assume that is a tightly bound s- 
wave state much smaller in size than the characteristic trap 
oscillation amplitude d,-: xif±, d,n) - uoir±, d„i) with d„i <K dr. 

However, the eigenstates of /if are not exact eigenstates 
of the total closed channel because they do not include the 
interaction of the spectator atom with the dimer in the state 
X- Nevertheless, in the regime of a tightly bound dimer, this 
residual atom-dimer interaction is expected to be small due to 
Pauli suppression. Suppose that atoms 1 and 3 have formed 
a dimer and atom 2 is the t spectator. The u„{s+,dr) part 
of the wavefunctions indicates that atoms 1 and 2 are sep- 
arated at a distance of the order of (3n)'^^t/r/2. Moreover, 
since they are fermions, Pauli exclusion means that their rel- 
ative wavefunction must have a node at rj = r^,, so that u„ 
remains small until they are separated by a length scale set 
by the trap length d,-- Since atom 3 stays close to atom 1 in 
a tight dimer, this means that atoms 2 and 3 also stay apart 
at a distance of the order of dr where the interaction V^' be- 
tween them is weak. Therefore, we can, to a first approxima- 
tion, neglect this spectator-dimer interaction so that approx- 
imate eigenstates for the closed-channel Hamiltonian /i"^' are 
antisymmetrised combinations of the two possible spectator- 
dimer states: 

[u„(s+,dr)x(r+,d„,)- u„(s-,d,.)xir-,d„,)] 

V2 

(9) 

for n = 0, 1, . . ., with energies (n + 1 12)tujj + £h + Sb- 

We can now discuss the full Hamiltonian with Feshbach 
scattering between the open and closed channels. The Hamil- 
tonian for the three atoms is given (in the centre-of-mass 
frame) by: 



= [lT,T,i> IT,T,°)] 



<T,T,il 
a,T,°l 



(10) 



where a{r) describes the coupling between the open and 
closed channels for two atoms at a distance of r. The eigen- 
states of this Hamiltonian ([10]) are now superpositions of the 
open and closed channel states: 



(11) 



where the first term represents the three atoms in the open 
channel in terms of the basis set Q and the second term rep- 
resents the atom-dimer states in terms of the basis set (|9|. 
An eigenstate with energy E obeys the Schrodinger equation: 
//TTi|(^TTi^ - /ili^TTi^^ This requires 

^ ^akinttki + b„{E„ + v) = Eb„, 



^ aunbn + {sk + si)au = Eau, 



(12) 



where e„ - {n + \ /2)hcL), ai.i„ - {kl\W\Xn) is the Feshbach scat- 
tering matrix element, and v - 6b + et is a detuning parameter 



which sweeps across the Feshbach resonance as a function of 
applied magnetic field. 

In the absence of the confinement potential, the resonance 
occurs when the atom-dimer state of energy + 6b coincides 
in energy with the zero-momentum scattering state (energy 
zero) in the open channel. This corresponds to v = 0. In the 
presence of the harmonic trap, the Feshbach scattering may 
cause resonances whenever the open and closed channel states 
are degenerate: 



{k + l+ \)nu) = V + (« -H 1 12)na) 



(13) 



for non-negative integers k, l{+ k) and n. A non-zero Feshbach 
coupling will lift the degeneracy and one will find a level anti- 
crossing in the energy spectrum as a function of v. We will see 
in section III that there are level crossings allowed by selection 



rules. 

To simplify our calculation, we can eliminate the open- 
channel coefficients a^i using the second equation in ( [T2| . 
Then, we obtain a matrix equation for just the atom-dimer 
components, bn, of the eigenstate: 

(E - E„)bn + M,„„(E)b„ = vb„, (14) 



where 



2 Ek + 



OklnO!klm 



Ek+Sr 



(15) 



The matrix M„,„ can be interpreted as representing an effec- 
tive scattering of the closed-channel state \xm) to l;^',,) via the 



virtual state \kl) in the open channel. We solve equation ( 14 1 



numerically by treating it, at any given E, as a matrix eigen- 
value problem for eigenvalue v. 

We need to evaluate the Feshbach matrix elements Ukin- We 
will assume that it only occurs when two atoms are close to- 
gether at a scale much smaller than the trap length dy We 
see that there are direct and exchange terms in this matrix el- 
ement. The direct term involves a(?"±)^(r±, d,„). For a tightly- 
bound dimer in the closed channel, we can reduce this to 
an effective delta-function coupling between the open and 
closed channelllSl [191 : a{r)x(r,d,n) — > a6{r). This leaves 
Of as a single parameter controlling the strength of the Fes- 
hbach interaction. There is also an exchange term involving 
a(r^)Unis-t, dr)x{r±, dm). This term is only large when all three 
atoms are very close to each other. These terms can be ne- 
glected due to Pauli suppression by the same argument that 
we neglected atom-dimer interactions in the closed channel: 
the range of the Feshbach scattering is much smaller that the 
typical separation (~ dr) of two atoms in the same hyperfine 
state. 

With this approximation of a short-ranged Feshbach inter- 
action, the matrix elements M„,„(E) can be written in terms of 
the Green's functions of the system: 

MnmiE) = —Idxdxdydy ij/„(x,y)ij/jx ,y ) (16) 
X [G^^\x,y;x',y';E) - G'^^\x,y;y',x';E)] 
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where ij/n - u„{s+,dr)5{r+) — Un{s-,dr)5{r_) and G*^' is the 
two-particle Green's function of non-interacting atoms in a 
ID harmonic well 



G^^\x,y;x',y';E) = ^ 



Ukix, d,)ui(y, dr)uk(x',dr)ui(y', dr) 
Sk+si- E 

(17) 

We notice lil2J that the harmonic oscillator Hamiltonian /!"^^p 
is invariant under rotations in the x-y plane and that the (r-t,5-t) 
are both related to (x,y) by such a rotation: r-t - [( VS + l)x -H 
(V3±l)y]/2V2andi± = [(1± V3)x-i-(l+ V3)y]/2 V2. There- 
fore, G^^\x,y; x'y'; E) = G'^^'(r-t, s±; r^, s'j_;E). This allows us 
to express all the integrations in terms of r and s variables. 
We find that the matrix elements M„,„{E) can be written in the 
final form 



M„JE) ^ a^i^G^'\0,0;E - s„)6„. 



(18) 



- IdsG 



-,0;E-s„ 



U„{^-^,dr^U„{s,dr)\, 



where G''^ is the single-particle Green's function 
G^^\x,x';E) — YinUn{^^dr)u„{x' ,dr)l{B„ — E). The ad- 
vantage of expressing the matrix in this form is that we know 
the quantities we need analytically |20|: 



G*"(x,0;£) = 



yfn d,- ha) 



r V2i 



dr 



.(0), 



(19) 

where Dy(x) is the parabolic cylinder function with Dvdyl) - 
2v/2gV/4{/(_y/2, i/2,3;2/2) where U is the confluent hyper- 
geometric function and DyiO) = 2"'^ ^/^/^ [(1 - v)/2]. 

The delta functions in the il/„tffm factors in the form ( [T6| in- 
dicate which particles are interacting via the Feshbach scatter- 
ing. The two G*^^' terms comes from the fact the virtual exci- 
tation via an open-channel state involves a direct process and 
an exchange process. The direct process gives the first term 
in ( [T8| which requires r+ - - ri - or = - r2 = 0: 
it involves only the interaction of one pair of atoms while the 
third atom remains a spectator. The exchange process involves 
the virtual dissociation of a dimer into t and X atoms in open- 
channel followed by the X atom recombining with a diff'erent 
t atom to form a new dimer. This gives the second term in 



( 18 I which is a term present in the three-body problem but not 
in a two-body problem. 

Before proceeding to discuss our results, we should point 
out the difference between | t, T> i) and I i, i, T) states. Since 
the Feshbach interaction only affects the X state in the open 
channel, there is no symmetry under the operation T*-*-l- The 
difference is that, in the XXT case, there is no residual in- 
teraction between the |o) atom in the dimer and the specta- 
tor X atom: the closed-channel attraction V^' only acts on a 
I t °) pair. However, we have already dropped this interac- 
tion within our approximation of a tightly bound dimer state. 
Therefore, our results should apply also to the XXT case. More 
explicitly, for the XXT case, the Hamiltonian for the relative 



degrees of freedom is of the form 



/z°P W+ W- 
W+ /i^ 
W- h"} 



V'-. 



(20) 



where W± = a(r±), V = [ | X,X,T) l°,X,T) I X,°,T) ] and 
and h^l have been defined in (j7]i. Now, the wavefunctions 

u„{s^,dr)x(r+,d,n)\°,i,'\) - Un{s-,dr)x{r-,d,y,)\ X,°,T) (21) 

are exact eigenstates of the closed-channel part of the Hamil- 
tonian. 

III. SELECTION RULES 

We can take advantage of selection rules for the Feshbach 
coupling. In the pseudospin language of T and X states of the 
open channel, we are discussing the three-atom states with 
S , - +1/2. If we completely antisymmetrise the wavefunc- 
tion of the three atoms EH : 

|«A(1, 2, 3)}a = |.A(1, 2, 3)) + 1^(2, 3, 1)) + \ip(3, 1, 2)>, (22) 

we can classify the three-atom states as S - 3/2 or 1/2. The 
S =3/2 spin wavefunction is totally symmetric under ex- 
change and so its spatial part is totally antisymmetric under 
exchange, leading to a node at - ri^2- Therefore, this state 
is unafifected by the .s-wave Feshbach scattering. On the other 
hand, the 5 = 1/2 states have mixed exchange symmetry (see 
for example [22 1) and Feshbach scattering is allowed. To see 
this explicitly in our formulation of the problem, we can write 
the open-channel part of the Hamiltonian in terms of the polar 
coordinates ^ = (x^ +y^y^^ and = arctan(y/x). In this coor- 
dinate system, the energies can be written as {2nr + mi + VftuD 
with integers > and m/ > 0. The 5=3/2 states corre- 
spond to mi - (mod 3) and these states have zero scattering 
matrix element to the atom-dimer states. The 5=3/2 state of 
lowest energy has energy Atio). 

In addition to exchange symmetry, we note that the system 
is symmetric under inversion {rxxi — * -'"i.i.s)- This means 
that the eigenstates of the system must be even or odd under 
inversion. In terms of the matrix equation ( [T4| , the solutions 
divide into two sectors, one for b„ with even n and one for 
b„ with odd n. More precisely, the matrix elements aun are 
non-zero only \fk + I and n are both even or both odd. 

There are also other matrix elements auin which vanish be- 
tween the open and closed channels: it can be shown that aun 
vanish when n > k + I. This can be proved using the fact that 
x'e^-'^^ H„{x)dx - for r < n ll24l . In other words, atom- 
dimer states with energy v-i-(n-i- 1 /2)ha> do not couple directly 
to the three-atom states with the relative energy (k + I + l)fiLo 
when n > k + I, even if they are degenerate (possible when 
V < hcj/2). However, they can be coupled by multiple Fes- 
hbach scattering. Therefore, we expect avoided level cross- 
ings for such levels but the splitting at the anti-crossing point 
will be smaller in size than the level splittings at anti-crossing 
points for levels coupled directly by a non-zero matrix ele- 
ment ffj./,,. This is because the splitting will be higher order in 
a. 
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IV. TWO-ATOM SYSTEM 

In this section, we will outline the results for the Feshbach 
resonances of two atoms in a ID harmonic well. This will be 
useful for the discussion of our results for three atoms. These 
results are analogous to the 3D results of Diener and Ho fSj. 
Our calculation differs from |8| in the lower dimensionality 
and also in the fact that we have factored out the motion of the 
centre of mass. 

In the centre-of-mass frame, the energy for two Feshbach 
interacting atoms in a well within the two-channel model is 
given by the equation 



^{l 2L) 



2hu)dy 



inul 



(23) 



The lowest-energy solution of this equation has an energy be- 
low fui>/2 for any value of the detuning parameter v. Since 
Ha)/2 is the lowest energy of the open-channel scattering 
states, this corresponds to a bound state. The existence of a 
bound state for any v is not unexpected in one dimension. 

We can relate our two-channel parameters a and v to the 
parameters of a one-channel scattering theory in the absence 
of the trap: aid, the i-wave scattering length, and rio, the 
effective range. These are the physical parameters that can be 
found experimentally. To do so, we follow the procedure of 
Diener and Ho JSl and examine the bound state energy of two 
atoms in the limit of an infinitely shallow well: w — > at a 
constant negative E. We also need \E\ v so that the state 
is clearly separated from dimer states (energy v) in the closed 



channel. In this limit, the equation (23 1 becomes 



yuor 
Tfiv 



lxa^\E\ 

ff.y2 



1/2 



0. 



(24) 



On the other hand, we find that the energy of a weakly bound 
state in a homogeneous case with no external confinement is 
given by the equation I/oid - A""id|£^I/^^ ~ (2a'|£^I/^^)'^^ = 
where \E\ - h^K^/{2fi) and /c/? <K 1 with R being the range of 
the scattering potential. 

Comparing those two equations, we find 



l/fliD = i^a^Kti^v) 



riD 



(25) 



We see that the effective ID scattering length is propor- 
tional to the detuning parameter v in the regime of a weakly- 
bound state where \E\ <K v. In terms of the parameters of the 
two-channel model, this regime corresponds to a detuning of 
V » (jja* /2h^y^^ . Note that this ID case is very different to 
the 3D case where 03/5 1/v, see section [Vl| 



V. RESULTS FOR THREE ATOMS 

We will now discuss our results for the spectrum and eigen- 
states of the three-body system. As mentioned already, we 



treat the matrix equation (14 1 as an eigenvalue problem for 
the detuning parameter v at fixed energies E. 

In order to obtain numerical results, we truncate the matrix 
equation ( [l4| by taking into account only atom-dimer energy 
levels with s„ < S2„c ■ This should properly reproduce the low- 
energy behaviour of the system for the few lowest atom-dimer 
energy levels. Here we choose nc = 20. We have checked that 
the ground and first excited states do not change significantly 
if we use a higher cutoff n^. 

We will discuss two opposite regimes of weak and strong 
Feshbach scattering. The kinetic energy of the system is on 
the energy scale ftoj. This should be compared to the magni- 
tude of the scattering matrix elements o-^fa. We note that the 
oscillator wavefunctions u„{x, dy) have a magnitude of the or- 

—1/2 1/2 

der of dy over a spatial size dy so that aki„ ~ a/dy . There- 
fore, a dimensionless measure for the strength of the Feshbach 
scattering is the ratio of these two energy scales: 



(26) 



tlOjdr 



1/2 



We will consider first the weak-scattering regime, a <«; 1, 
where the open and closed channel states are weakly coupled. 
Then, we will discuss strong scattering, a » 1, relevant to 
^"K atoms |[T8l|23l. 

We can also obtain a scale for the range of detuning 6v over 
which the Feshbach scattering significantly affects the atom- 
dimer energy levels, as governed by (fT4ji. This can be es- 

' ' 1/9 9 

timated by the scale of the elements of M„„, ~ (a/dy ) /tiLO. 
Therefore, a dimensionless measure of the detuning parameter 
is vhaidy/a^. From the results (25 1 for the two-body problem, 
we see that this is also the ratio of the effective ID scattering 
length, fliD, for states in the open channel (in the absence of a 
trap) compared to the trap oscillation length dy: 



fllD 

dy 



vfuijdy 



(27) 



in the regime of a weak two-body bound state. 

We present in figure [T| our results for the spectrum at a = 
0.5. We plot the energy in units of Hai versus the dimension- 
less detuning parameter vha>dy/a^. The selection rules dis- 
cussed above are reflected in the energy spectrum as it can be 
seen in figure[T] The two graphs, (a) and (b), correspond to the 
energy levels for states with even and odd inversion symmetry. 

The 5=3/2 open-channel states which are unaffected by 

the Feshbach scattering are found at energies E/Hll) = 4, 6, 

At negative detuning and at negative energies, we see the un- 
perturbed atom-dimer states the different diagonal lines 
corresponding to different atom-dimer oscillations. As we in- 
crease V, Feshbach interaction with open-channel states be- 
come possible near the resonance points given by ( [T3| ). At the 
resonances, there are anti-crossings between the three-atom 
and atom-dimer states. For small a, we can approximate by 
restricting the matrix equation ([14]) to just the two open and 
closed levels near the anti-crossing in question. Then, we see 
immediately that, for resonances at positive detuning v, the 
size of the splitting at resonance is of the order of the matrix 

1 /2 

element coupling the two states: atin ~ a/dy - aha>. 
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However, for resonances at negative v, the splittings are 
much smaller. As discussed in the section [III] these are the 
levels which are only coupled through multiple Feshbach scat- 
tering. 



flOJ 
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FIG. 1: Energy levels (in the centre-of-mass frame) versus detuning 
for three atoms in a harmonic well with .s-wave Feshbach scattering 
at a dimensionless coupling a = 0.5. Dashed lines represent total 
pseudospin 5=3/2 states which are unaffected by i'-wave Feshbach 
scattering. Figures (a) and (b) represent eigenstates with odd and 
even inversion symmetry respectively. Grey circles represent eigen- 
states in the absence of Feshbach scattering. Black dots represent 
5 = 1/2 states. 

We can now turn to the regime of strong Feshbach scat- 
tering. Large coupling strength a means much stronger level 
anticrossing than for a - 0.5 case shown in figure [T] As an 
example, we have calculated the spectrum for a - 100. (We 
will see later that there is very little admixture of high-energy 
level in the lowest-energy eigenstates, justifying our use of a 
relatively low cutoff tic-) The six levels with the lowest ener- 
gies are shown in figurej2] The apparent kink in the highest 
branch in figures |2] and [3ps due to the level crossing that can 
be better seen for the small a case in figure [T] Again, we see 
the unperturbed open-channel states at ATitij and 6tLco (dashed 
lines). Note that the lines corresponding to the unperturbed 
atom-dimer states would be very steep and very close to the 
vertical axis on this graph. 

We find level crossings in the low-energy spectrum between 
these two sets of states with opposite parity under inversion. 
This is most easily seen in energy spectrum at weak Feshbach 
scattering (Fig.[T]). For large positive detuning, the lowest en- 
ergy states belong to the odd-parity branch that gives the open- 
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FIG. 2: Lowest six energy levels at strong Feshbach coupling: a = 
100. Dashed lines represent total pseudospin 5 = 3/2 states. Grey 
and black dots represent 5 = 1/2 states with odd and even inversion 
symmetry respectively. 



channel |01) state asymptotically far away from the resonance 
(v +oo). On the other hand, for large negative detuning, the 
states of lowest energy belong to the even-parity branch that 
becomes the \xq} atom-dimer state asymptotically (v — » -oo). 
Therefore, the two parity branches must cross as a function 
of detuning. For weak Feshbach scattering, the crossing for 
the lowest two parity states occurs at a positive detuning, just 
below E = 2tia>. For strong Feshbach scattering, we find that 
these level crossings have shifted to small negative detuning 

- 7 1/2 

and negative energies (vtitjdr/a ^ -0.07 and E ^ -9a/di. 
for a - 100). Similar features were found for the spectrum 
of three atoms in a 3D harmonic well within the one-channel 
model lITIl . 

The level crossings between the branches of opposite in- 
version symmetry are best seen if we subtract from the three- 
body energies (Ej) both the ground-state energy of two atoms 



in a well (E2 as given by (23 1) and the energy of one atom in a 
separate well {hull). This quantity, E^- E2-fujL>l2, is plotted 
in figure [3] The subtraction does not alter the crossings of the 
three-body states at any given v. In section VII we speculate 



that these crossings may lead to quantum phase transitions. 

This data allows us to discuss the ground state of the system 
with three atoms for every two wells. One might ask whether 
the Feshbach interaction could generate an attraction such that 
there is a regime where the system would prefer to arrange 
three atoms in one well and none in the other well. Such a 
distribution would have an energy of £3 plus Ticojl for the 
zero-point motion of the centre of mass. We can also have 
two atoms in one well and one atom in the ground state in the 
other well. This has energy £2 plus Sw/2 for the centre of 
mass in one well and SdLt/2 in the other well. Therefore, the 
energy gap between these two states is E^ - E2- fia>/2. From 
Fig. [3] we see that this gap remains positive for all detunings, 
suggesting that the ground state of an optical lattice with 3/2 
atoms per site to be a phase with uniform density. This energy 
gap has a minimum as small as Q.lhoj at a = 100. This mini- 
mum is due to a reduction in £3 compared to E2 arising from 
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FIG. 3: Difference between the energy of three atoms in a harmonic 
well and the ground-state energy of two atoms in the well plus the 
lowest one-atom energy in a separate well. Grey and black dots rep- 
resent 5 = 1/2 states with odd and even inversion symmetry respec- 
tively. 



exchange processes which do not occur for two atoms, as dis- 
cussed in previous section after equation ( [T9| . We note that 
this energy gap is much smaller in magnitude than the Fesh- 

1 /2 

bach energy scale a/d/ and the individual energies £3 and 
El- 

We can similarly ask about the ground state of four atoms 
in two wells. The evenly-distributed state with two atoms per 
well is lower in energy because it contains two binding ener- 
gies 2E2 whereas the uneven distribution can only take advan- 
tage of Feshbach physics for one pair of atoms in £3. 

Having discussed the energy spectrum of the three-body 
system, we will now examine the wavefunctions of the low- 
energy eigenstates. The compositions of the ground state and 
the first excited state (for pseudospin 5=1/2) are shown at fig- 
ure]?] Figure]4](a)/(b) presents the composition of the lowest 
energy state with an odd/even inversion symmetry. As ex- 
pected, the state in figure (a) becomes the atom-dimer 
level at large negative detuning while it is predominately the 
open-channel |0, 1) level at large positive detuning. Near the 
resonance, there is also a significant admixture of other open- 
channel levels and a smaller admixture of other atom-dimer 
levels. The most significant components are shown at the fig- 
ure. Figure (b) gives an example with even inversion symme- 
try. The state is asymptotically the \xo) atom-dimer state at 
large negative detunings. Closer to the resonance at negative 
V there is also a significant Ixi) component. For positive v, 
the dominant contributions are from the |0, 2) and 1 1, 3) open- 
channel levels with the |1,3) contribution vanishing far from 
the resonance. 

We notice that, while only \xi ) dominates in the lowest en- 
ergy state with odd inversion symmetry, both \xo) and Ixi) 
are significant for the corresponding state with even inversion 
symmetry particularly near the resonance. This difference can 
be traced back to the Feshbach matrix element cKi/„. As dis- 
cussed in]III] the wavefunctions involved should have similar 
symmetries to ensure a finite value for aki„. Thus, using the 
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FIG. 4: Ground and first excited pseudospin 5 = 1/2 state compo- 
sition as a function of detuning v (a = 100). The probabilities la^/p 
and \b„\^ correspond to open-channel and atom-dimer probabilities 
respectively (see the equation l]TTJ for definition of a^i and b„). (a) 
Lowest-energy state with odd inversion symmetry (lowest gray data 
in figure ]2j. (b) Lowest-energy state with even inversion symmetry 
(lowest black data in figure]2j. 



rules that k I (Pauli principle) and that k + I and n should 
be both even or both odd with n < k + I, the lowest few aki„ 
that contribute to the lowest energy states are q-qii and aoB 
for the odd-parity state and ao20 and 0-022 for the even-parity 
state. These matrix elements are large only if the quantum 
numbers k and I are similar to n and so: laonl » laonl and 
|Q'02oI ~ |Q'022I- In terms of equation ( 14 1 for the coupling of 
the closed-channel amplitudes, the magnitude of Mu is big 
compared to that of Mi 3 for the odd case, while both Mqo and 
M()2 are significant for the even case. Hence, both \xo} and 
IXi) contribute to the lowest energy even sector eigenstate in 

Fig.]4];b). 



These results allow us to understand the validity of our trun- 
cation of the matrix equation ( ]T6| . We see that the admixture 
of high-energy levels is small. This allows us to ignore the 
contribution of \x„) levels with n > 2nc = 40. For compar- 
ison, if we use a smaller cutoff of nc = 5 instead of 20, the 
value of V at a given E changes by less than 1 % and the ratios 
between different amplitudes change by less than 10%. 
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VI. COMPARISON WITH ONE-CHANNEL MODEL 

In this paper, we have used a two-channel model for states 
near the Feshbach resonance. It is often convenient to use a 
one-channel model in which the Feshbach-induced interaction 
between open-channel states is modelled by a contact poten- 
tial: V{r) - g6{r), with a potential strength g that depends 
on the de tuning liTj J2 1 . We discuss now how a one-channel 
model may be appropriate for the three-body problem in a ID 
well. 

Physically, we expect that an effective model for just the 
open channel cannot capture physics involving the details of 
the dimer state in the closed channel. Therefore, one-channel 
models can approximate the results of two-channel models 
only if we focus our attention on scattering states and weakly 
bound states, neither of which would have a strong admix- 
ture of the closed-channel dimer. We have been discussing 
the branch of bound states induced by Feshbach interaction. 
In three dimensions without a harmonic trap, these states are 
weakly bound, with a binding energy that vansihes as we ap- 
proach the resonance (v - 0) from negative detuning. In other 
words, these eigenstates have only a small amplitude in the 
closed channel dimer state. Therefore, the matching to a one- 
channel model is best near the resonance. The range of va- 
lidity of the one-channel model will be wide if the Feshbach 
coupling is strong. 

This comparison between one-channel and two-channel 
models in 3D has also been studied in the presence of a har- 
monic trap. At a wide Feshbach resonance (strong coupling 
between the open and closed channels), the admixture of the 
closed channel dimer component to the low energy states is 
very small 1 8 1 and the effective interaction between the atoms 
close to the resonance can be described by an effective contact 
potential ||2l[l2]|25]|26l withg3/j = 2;r?!^a3D/yU which diverges 
as one sweeps the system across the resonance at v = 0. How- 
ever, there are also narrow resonances where a one-channel 
model does not apply |27 1 due to the significant admixture of 
the closed-channel component. 

We will now turn to our results in one dimension. We have 



akeady studied in section IV the matching of the two-channel 



results (for two atoms without a trap) to one-channel scatter- 
ing parameters for the branch of bound states that develops at 
energies below the open channel. Note that this branch is not 
confined to v < (as is the case in 3D without a trap). In- 
stead, it exists for all v with the most weakly bound states at 
large positive detuning. This is the regime where we expect 
the one-channel and two-channel models to converge. The 
weak-binding condition that l^l ~ tf/2fia^^ <g: v can be ex- 
pressed in terms of the dimensionless Feshbach and detuning 
parameters used in our study ( [25| ): 



vftCOdr a ID 1 



(28) 



We see that this matching applies to positive detuning away 
from the level anticrossing near v = 0. Note that the regime 
of validity improves at strong coupling: a ^ oo so that such a 
theory may be applicable for realistic '^''K systems except very 
close to the resonance. In terms of an effective contact poten- 



tial strength, we have giu = -fP'/i^air, fl] so that the regime 
of validity corresponds to a weak contact potential which traps 
a bound state for arbitrary weak gi^ in one dimension. 

These conclusions about the validity of a one-channel ap- 
proximation are also consistent with our three-atom results 
presented above. We can examine directly the degree of the 
admixture of open- and closed-channel components as a func- 
tion of detuning. At large negative detuning, the states illus- 
trated in figure [4] are predominantly the closed-channel atom- 
dimer state and so should have no relation to any model based 
on open-channel states. The change in the balance between 
open- and closed-channel components as a function of in- 
creasing detuning is clearly seen in figure [4] Indeed, in our 
ID system, there is a significant admixture of the atom-dimer 
levels close to the resonance even for strong coupling between 
the open and closed channels {a » 1). Therefore, we con- 
clude that the system cannot be modelled by a one-channel 
theory close to the resonance. On the other hand, for large 
positive detuning, the states are predominantly open-channel 
states and so can be modelled in a effective theory for the 
open-channel states. 




0.25 



FIG. 5: The lowest energy state with odd inversion symmetry for 
three atoms calculated using the two-channel models (solid circles) 
and one-channel model (open circles). The horizontal axis is plot- 
ted using the correspondence in equation l|28j. Grey line represents 
the closed-channel state in the two-channel model in the absence of 
Feshbach coupling (E = v + 3ha)/2). 

Furthermore, we have performed a one-channel calculation 
for three atoms in a well for our system using an effective 
^iD 128], following an analogous study for the 3D problem 
lfT2ll . From Fig. |5] we see that the energies calculated for 
the lowest branch of eigenstates agree in the one-channel and 
two-channel models at large positive detuning, but they differ 
significantly close to the resonance. Similar conclusions can 
be made for the two atoms in a ID well by comparing one- 
channel results f71| the two-channel result ( [23| l |8|. We can 
also confirm that the one- and two-channel calculations give 
different results when there is a significant admixture of the 
atom-dimer levels. 

Interestingly, we find that the energy difference £3 - £2 are 
very similar in the one-channel and two-channel models for 
the whole range of positive detuning at a = 100. 
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VII. SUMMARY AND DISCUSSION 

We have considered two species of feimionic '^''K atoms 
at ultra cold temperatures trapped in a ID optical lattice and 
close to an inter-species Feshbach resonance. We have con- 
centrated on the zero-tunnelling limit in which the dynamics 
of the system is well described by the local spectrum of atoms 
at a lattice site. The calculations are performed within the 
two-channel model for two and three atoms at a lattice site 
approximated by a harmonic well. 

Close to the resonance, we find a significant admixture of 
the atom-dimer states to the lowest energy eigenstates of the 
system. Consequently, the one- and two- channel models must 
differ close to the resonance for our ID system since a one- 
channel theory cannot capture the physics of the dimer state. 

The results suggest that the ground state is a state with uni- 
form atom density. In particular, for three atoms in two wells, 
we compared the energy of all three atoms in one well (3-1-0) 
to the energy of two atoms in one well and one atom in the 
other (2+1). The energy difference can be as small as Q.ltuD 
at a = 100 which is much smaller than the Feshbach energy 
scale. However, we note that, for the S, - +1/2 case, we 
have neglected the residual interaction between the atom and 
dimer in a [;^'„) state. There is an attractive component from the 
closed-channel potential V^' and also a Feshbach scattering 
term, both of which we argued to be suppressed by fermionic 
statistics. We can ask whether these additional terms may sig- 
nificantly reduce or even change the sign of the energy dif- 
ference between 2-1-1 and 3-1-0 configurations. In other words, 
we ask whether there is an instability to a density wave in an 
optical lattice with a filling close to 3/2 atoms per site. 

To obtain a quantitative analysis for such a scenario, we 
would require a detailed short-distance description of the 
dimer state. Here, we will limit ourselves to a rough estimate 
for the lowest-energy state at positive detuning by considering 
the probability that the spectator atom can be in the vicinity of 
the dimer 

This state has odd inversion symmetry and the amplitude 
for the two atoms (r^ and r2) to approach each other is pro- 



portional to (ri - r2)/d/ . So, the probability for the spectator 
to be in the vicinity of a dimer of size dm is proportional to 

j "'ir/d,-)^dr/dr ~ (d„,/dr)^. We estimate the residual attrac- 
tion fad between the atom and the dimer would be of the order 
of V{d„,ldrf' where V is a measure of the strength V"^' expe- 
rienced by the closed-channel dimer, e.g. the minimum value 
of the potential. As mentioned above, the minimum energy 
gap in £3 - £2 - tL<x>l2 is only a fraction of tio). Therefore, 
if V is sufficiently large compared to ticj, there is a possibility 
that this additional attraction can favour the 3-1-0 configuration 
compared to the 2-1-1 configuration in two wells. (A similar 
term in the Feshbach scattering can be absorbed into a renor- 
malised scattering strength a, but our results are not sensitive 
to the value of a in this regime of large a.) In fact, the even- 
symmetry state has a similar but larger correction of Vdm/dr 
which could change the ordering in energy of the even and 
odd branches over a range of detuning. These speculations 
need to be confirmed by further work. 

Note that these hypothetical scenarios rely on the residual 
interaction between an t atom (but not a i atom) with a dimer 
in the state | t °)- In other words, such a scenario would 
reveal the asymmetry of the system under the transformation 
t<->J,, it would only be a possibility if we have more atoms in 
the t open-channel state than in the i state. The magnitude 
of this effect depends on the ratio of the size of the dimer 
in the closed channel compared to the trap size and therefore 
depends on the details of the actual Feshbach resonance. 

Finally, we discuss briefly the optical lattice with weak tun- 
nelling between a set of wells. Consider the case with three 
atoms per site on average. From the level crossing we found 
in figure [3] we see that the ground state would consist of states 
with odd or even inversion symmetry within each site depend- 
ing on the detuning parameter. As we sweep past this cross- 
ing, this may be a first-order transition. However, it is also 
possible that a continuous quantum phase transition would oc- 
cur because an inter-site tunnelling would couple the states of 
different inversion symmetries. We need to investigate collec- 
tive quantum fluctuations to study this possibility. 
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